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We investigate several quantifiers of the electroencephalo- 
gram (EEG) signal with respect to their ability to indi- 
cate depth of anesthesia. For 17 patients anesthetized with 
Sevoflurane, three established measures (two spectral and one 
based on the bispectrum), as well as a phase space based 
nonlinear correlation index were computed from consecutive 
EEG epochs. In absence of an independent way to determine 
anesthesia depth, the standard was derived from measured 
blood plasma concentrations of the anesthetic via a pharma- 
cokinetic/pharmacodynamic model for the estimated effective 
brain concentration of Sevoflurane. In most patients, the 
highest correlation is observed for the nonlinear correlation 
index D* . In contrast to spectral measures, D* is found to 
decrease monotonically with increasing (estimated) depth of 
anesthesia, even when a "burst-suppression" pattern occurs in 
the EEG. The findings show the potential for applications of 
concepts derived from the theory of nonlinear dynamics, even 
if little can be assumed about the process under investigation. 



I. INTRODUCTION 

Several recent publications have used analysis methods 
derived from the theory of nonlinear dynamics to describe 
aspects of neuronal activity. Measures like the estimated 
largest Lyapunov exponent (LLE) or the correlation di- 
mension (D2) were repeatedly applied J|J^| to quantify 
brain electrical activity, for example under the influence 
of epileptic processes, during different mental states, or 
in different sleep stages. Attempts to interpret the abso- 
lute numbers obtained in terms of "system complexity" , 
"attractor dimension" , "chaoticity" , or the like have met 
harsh and justified criticism ||. Any such interpreta- 
tion would require some narrow assumptions (station- 
arity, predominant low dimensional determinism, etc.) 
which are expected to be violated when biological time 
series are analyzed. On the other hand, commonly ap- 
plied linear or power spectral measures are based on dif- 
ferent, but equally narrow assumptions, including (weak) 
stationarity and Gaussianity. The success of these latter 
methods is undisputed because of the practical usefulness 
of the results. In this paper, we will try to take such a 
pragmatic attitude and compare a number of indicators 
of human electroencephalogram (EEG) time series with 
respect to their power in monitoring depth of anesthesia. 



Keeping patients at a well defined level of anesthesia 
is still a difficult problem in clinical practice. If anesthe- 
sia is too deep, a decompensation of the cardiovascular 
system is threatening. When anesthesia is too flat, the 
patient may wake up. Depth of anesthesia is expected to 
be reflected in the EEG. Some of the underlying mecha- 
nisms are known from animal studies. The EEG is gen- 
erated by electrical discharges (i.e. postsynaptic poten- 
tials) of neurons located near to the surface of the brain 
(predominantly neurons from the neocortex) H. Dur- 
ing the awake state, complex tasks can be processed by 
these cells. Neuronal activity changes during anesthe- 
sia as well as during "synchronized" sleep (including all 
sleep phases except for the rapid-eye-movement-sleep, in 
which most dreams are experienced). In the middle of 
the brain (i.e. in the thalamus), neuronal populations 
start to produce an oscillating activity that increasingly 
synchronizes neo-cortical neurons. Simultaneously, the 
excitability of neo-cortical neurons decreases so that the 
amount of interactions between them decreases as well. 
Finally, nearly all neo-cortical neurons fall silent. This 
state will not be reached during normal sleep, it is only 
observed during coma or deep anesthesia. A small num- 
ber of thalamic oscillators can even be active during this 
phase. Before a continuous zero line occurs in the EEG, 
a so-called burst-suppression phase can appear, where 
volleys from these thalamic neurons account for a cyclic 
excitation of the silenced neocortex ^ . 

The described features, except for the zero line, can 
be seen in Fig. |l|. Note the marked hysteresis: The burst 
suppression disappears only about 1 min after the plasma 
concentration of the anesthetic is reduced. The scenario 
may be summarized by the hypothesis that with increas- 
ing depth of anesthesia, the highly complex activities of 
the awake state are increasingly superimposed and re- 
placed by more simple oscillators. Of course, such a pic- 
ture amounts to a gross simplification. The full mech- 
anism will certainly not be well described by a single 
parameter like "depth of anesthesia" . 



II. VALIDATION OF MEASURES FOR DEPTH 
OF ANESTHESIA 

In current clinical practice, one or a few channels of 
EEG are routinely displayed during difficult anesthesias. 
Since the attending personnel has to monitor several crit- 
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FIG. 1. EEG from a patient during anesthesia. The con- 
centration of the anesthetic drug was increased until burst 
suppression first occurred (here at about 5 min) and then 
reduced again. The EEG is given in uncalibrated A/D units. 

ical parameters (blood pressure, heart rate, etc.), the 
vast amount of information contained in the EEG must 
be severely condensed in order to be useful. Only a few 
numbers may be monitored at a typical intervention time 
scale. Most pragmatically, a single number should be 
produced that indicates the instantaneous depth of anes- 
thesia of the patient. 

In order to derive and to validate such a measure of 
depth of anesthesia, the "true" level of anesthetic depth 
has to be known. One possible definition is to take the 
concentration of the anesthetic drug in the brain tis- 
sue as an indicator. This is reasonable if a monotonic 
relationship between the concentration of a single drug 
and the depth of anesthesia can be assumed. However, 
when multiple drugs are used, the drug interactions may 
become unpredictable. Unfortunately, even if a single 
drug is given, its brain concentrations cannot be mea- 
sured directly in human patients. In order to circumvent 
this difficulty, in this study the gaseous (volatile) drug 
Sevoflurane was used which can be employed as single 
anesthetic. In that case the concentration of the drug in 
the expired air can be measured, which is proportional 
to the blood-plasma concentration of Sevoflurane. Then, 
by considering the brain as a hypothetical effect com- 
partment, the equilibration of plasma concentration C p 
and the concentration in the brain C e can be modelled 
in a standard way by 

^ = k(C p -C e ) (1) 

where k is the constant of equilibration that will vary 
from patient to patient. Since k is unknown it has to 
be obtained by a fit of this pharmacokinetic model to 
the data. Equation (jl]) expresses the hysteretic response 
noted above to a change in the amount of Sevoflurane 
administered. 



As a further complication, the pharmacodynamic rela- 
tion between the concentration in the brain as the effect 
compartment C e and the neuronal effect E is not known 
a priori. A reasonable assumption is that there exists a 
one-to-one, monotonic relationship. In that case, and if 
C e (t) and E(t) are sampled at discrete times sufficiently 
far apart to assume statistical independence, this can be 
expressed by Spearman's rank correlation. Independence 
is not a reasonable assumption for consecutive values of a 
slowly varying quantity. Therefore we cannot use the ab- 
solute value of Spearman's rank correlation for a formal 
statistical test. However, since the conditions are identi- 
cal for all the methods we want to compare, significant 
differences in the correlation may be evaluated. 

The aim of our analysis here is to relate the measured 
C p and some estimated value E of the neuronal effect 
E, assuming the kinetic model, Eq. ([!]). Let Ri be the 
rank of C e (i) among all available C e 's and Si the rank of 
E(i) among all the corresponding i?'s. An estimate for k 
can then be found by minimizing the sum of the squared 
differences of ranks: 

N 

E(^-^) 2 - ( 2 ) 

i=l 

The corresponding rank correlation coefficient || 

= Ei(fr-i?)(ft-5) 

quantifies the success of the fit and can thus be used as a 
performance index for a method to estimate the neuronal 
effect. 

For our study, 17 patients undergoing an elective 
surgery were anesthetized with Sevoflurane. As described 
above, being unable to directly measure the brain con- 
centration C e , we measured the concentration in the 
breathed air at the end of expiration, which is propor- 
tional to the blood-plasma concentration C p . The neu- 
ronal effect E is estimated for each patient by four dif- 
ferent methods based on bipolar EEG recordings from 
frontal scalp regions. All EEG recordings used in this 
study were sampled at 256 Hz with 12 bit resolution. 
For the time resolved evaluation, a window length of 30 s 
was used with windows taken every 5 s, that is, with 83% 
overlap. The total observation period for each patient 
varied between 13 min and 55 min, median 34 min. 

III. EEG BASED MEASURES OF ANESTHESIA 
DEPTH 

In this study, we have compared four methods to derive 
an indicator for depth of anesthesia from EEG data. Two 
of these measures are based on the power spectrum, one 
additionally takes the bispectrum into account, combined 
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with a special feature detection scheme. The fourth is a 
nonlinear correlation index based on phase space recon- 
struction and correlation sums. Like in any other prob- 
lem of EEG characterisation, none of the used methods 
is derived from first principles nor can any one be fully 
justified on physiological grounds. Rather, they aim at a 
useful reduction of the vast amount of information in an 
EEG trace to one or a few numbers that may be moni- 
tored with justifiable effort. 



A. Power spectral measures 

Apart from visual inspection, the most common 
method to analyse and quantify time series data in clin- 
ical medicine — often the only viable one with the given 
data quality — is the frequency power spectrum, usually 
computed using the Fourier transform. Spectral estima- 
tion in itself provides only a moderate reduction of in- 
formation, except in the case that dominant oscillators 
can be identified by sharp spectral lines and overtones. 
The art of spectral analysis is to define key parameters 
that turn a power spectrum into a number, or at most 
a few essential quantities. For line spectra, the positions 
and widths of the resonances are obvious candidates. For 
broad band spectra like the EEG, one may consider the 
spectral content in certain frequency bands, characteris- 
tic frequencies, etc. 

Here we use two "classical" spectral parameters |?] || 
that have been previously applied for this purpose. The 
spectral edge frequency 95% |l0| indicates the maximal 
relevant frequencies while the median frequency [ p"T[ gives 
a rough indicator over the overall typical frequency of the 
recording. The physiological motivation lies in the expec- 
tation that at deeper levels of anesthesia, high frequency 
components of the surface EEG should be suppressed due 
to the entrainment by slow oscillations in more central re- 
gions of the brain. Therefore, both indices are expected 
to decrease with increasing depth of anesthesia. 

Both measures are based on the time windowed Fourier 
transform (using 5% cosine tapering of window edges) 
H U) = Et i Io 1/l ( < )exp(2 7 ri<//iV) of the digitally sam- 
pled EEG time series hit). Inspecting the spectral band 
between /i ow = 0.25 Hz and /high = 30 Hz, the accumu- 
lated periodogram 



(4) 



is used to determine the two spectral parameters, the 
spectral edge frequency 95%, here denoted by / se g5, and 
the median frequency /median- The parameter / se 95 is de- 
fined as the frequency at which P(/ SC 95) = 0.95P(/hi g h) 
and the median frequency / mc dian as the frequency at 

which P(/ m edian) = 0.5P(/ h igh)- 



B. Bispectral index 

A popular parameter in EEG instrumentation for anes- 
thesia is represented by the bispectral index (bis) [ [l2|| . 
This composite parameter uses multiple measures of the 
EEG power spectrum and bispectrum as well as burst- 
suppression analysis and has been optimized to predict 
the degree of sedation. Unlike all other EEG param- 
eters, it has been validated in large studies as a mea- 
sure of sedation in anesthesia 0]. Unfortunately, the 
precise algorithm is proprietary and has not yet been 
published. The bispectrum is one possible natural gen- 
eralisation of the Fourier power spectrum to nonlin- 
ear signals. While the ordinary spectrum may be de- 
fined as the Fourier transform of the two-point auto- 
covariance function (x(t)x(t — r)) of a signal x(t), the 
bispectrum is obtained by a two-dimensional Fourier 
transform of the three point auto-covariance function 
(x(t) x(t — Ti)x(t — T2)). The way the bispectral index 
condenses this two-variate function into a single number 
seems to have been established heuristically. At least 
no physiological motivation is given in the literature. In 
the absence of a published definition, we use the values 
of the bispectral index as computed by the Aspect 1000 
EEG monitor and record it every 5 s. 

The reported relative success 13 of the bispectral 
index suggests that properties of the EEG beyond the 
power spectrum are important for distinguishing of dif- 
ferent states of anesthesia. We do not want to enter into 
a discussion here of what nature the extra structure is. 
It is very difficult to disentangle the effects of nonlin- 
earity and non-stationarity, and, if nonlinearity is found, 
to what extent it is of dynamical nature or if it merely 
reflects the complicated transduction properties of the 
intervening tissue. From a physiological point of view, 
the situation is rather involved. The individual neurons 
clearly show nonlinear behaviour in their action poten- 
tials, but the law of large numbers suggests that this 
structure might be averaged out in surface EEG traces. 
Non-stationarity is very prominent in anesthesia, in par- 
ticular in the burst suppression phase, where episodes of 
low amplitude activity are interrupted by short outbreaks 
of large amplitude waves. Spectral indicators are known 
to give spurious values in the presence of burst suppres- 
sion. Therefore, the bispectral index algorithm contains 
a pattern recognition scheme to detect burst suppression, 
which is then treated separately. 



C. Nonlinear correlation index D* 

Early claims of low dimensional strange attractors in 
brain signals have been identified as spurious due to the 
lack of certain precautions when using nonlinear analysis 
tools. (See [| 14 for further discussion and references.) 
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This does not exclude that tools and methodologies de- 
rived in the context of dynamical systems, or chaos the- 
ory, may be still be useful when employed with care. In 
particular, interpretations in terms of fractality, chaotic- 
ity, complexity, and the like, have to be attempted with 
utmost discretion. 

In that spirit, we have adapted a prescription for 
an overall index of (nonlinear) coherence that has been 
found powerful for anticipating epileptic seizures from 
implanted electrode recordings |1| as well as in epilepsy 
models on a cellular level [ fl6| |. This index, which we 
will call D* contains many ingredients familiar from the 
Grassberger-Procaccia algorithm for the correlation di- 
mension Jl7],[l8|], whence the D. The * is there to remind 
us that there are important differences that preclude an 
interpretation in terms of fractal dimensions, number of 
degrees of freedom, etc. 

Starting form a scalar time series s ni n = l,...,N, 
we first form delay embedding vectors as usual: s„ — 
(s„_( m _i) T , . . . , s„_ r , s n ). The correlation sum C m (r) is 
then defined as usual by 



JV i-An 



C m (r) = aJ2 E 9 (r- IK - s iH) 



(5) 



where O is the step function and a = 2/ (N — An)(N — 
An — 1) is a normalisation constant. For the estimation 
of the correlation dimension, one would take An large 
to exclude temporal correlation and study the slope of a 
double logarithmic plot of C m (r) versus r and look for a 
scaling region, that is, a range of values of r where 



C' m (r) = 



dlog(C m (r)) 
dlog(r) 



(6) 



is constant. For high enough embedding dimension m, 
that constant would be an estimator of the correlation 
dimension Z?2- There is ample literature on using the 
correlation integral for dimension estimation ]3 . 

Except for specific pathologies, we do not expect this 
procedure to arrive at a reasonable, finite value for D 2 
when applied to the EEC Since we cannot hope to es- 
timate a proper dimension in the first place, there is no 
particular reason to concentrate on geometrical correla- 
tions only, or to aim at a proper phase space embedding. 
In principle, we could instead optimize parameters for 
the specific purpose at hand. However, a total of 17 pa- 
tients is too small a population to be split into a training 
set and a test set and we are not interested in in-sample 
results that may be due to over-fitting. Therefore, we 
copy the parameter settings and operational procedures 
from previous studies on epilepsy prediction ]T5| without 
further optimization. Using unit delay, a range of em- 
bedding dimensions m = 10, . . . , 25 was studied and only 
immediate temporal neighbours with i — j < An = 10 
were excluded. Each segment of 30 s was digitally band- 
pass filtered (0.25 Hz-30 Hz) before computing C m (r). 
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TABLE I. 


Results of the fit for the relaxation parameter k. 



The operational procedure of determining a "plateau" 
value D* of C' m (r) is also described in Ref. ]Tq] . It is 
worth noting that this automatic procedure is designed 
to exclude both, large amplitude artefacts (large r) and 
small amplitude fluctuations (small r). Some conjecture 
on the reasons for the success of this procedure will be 
offered in the discussion section below. 



IV. CLINICAL RESULTS 

In all 17 patients the following protocol was carried 
out. The concentration of Sevoflurane was at first slowly 
increased. As soon as a burst suppression pattern or 
even a zero line occurred in the EEG, Sevoflurane con- 
centration was decreased again. Two typical examples of 
corresponding EEG time series and estimated measures 
of anesthesia depth are shown in Fig. 0. In case B, im- 
mediately a zero line occurs in the EEG rather than a 
burst suppression pattern. It is important to remember 
that the actual neuronal anesthetic effect is not simply 
proportional to the plasma concentration C p shown in 
Fig. H but follows hysteretically through the relaxation 
phenomenon we modelled by Eq. ([j]). Generally speak- 
ing, all four measures studied are expected to show a 
negative correlation with the depth of anesthesia. 

In both cases, a slight decline of D*, but a small in- 
crease of bis precedes the burst suppression (resp. zero 
line) phase. The reaction of bis on these abrupt changes 
of EEG patterns is delayed. Depending on frequency 
characteristics of the background noise, spectral mea- 
sures either decrease (A) or increase (B) during phases 
with zero-line EEG. For a quantitative evaluation, we 
need to estimate the effective brain concentration C e 
of Sevoflurane. For this purpose, we assume the phar- 
macokinetic model given by Eq. (Q) and determine an 
equilibration-constant k individually for each case by 
minimising the sum of the squared differences of ranks 
between C e and E, c.f. Eq. (||). The values for k that 
were found [jl9| are summarized in Table |. 

Compared to the variation from patient to patient, 
there were no dramatic differences among the results for 
the different measures of anesthesia. It is important to 
notice that the fit for k is only important in order to es- 
tablish a monotonic relationship between an EEG mea- 
sure and assumed depth of anesthesia. The parameter 
k describes the time course of the equilibration of anes- 
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FIG. 3. Signed ranks of the differences between the cor- 
relations for D* and bis. A bar to the right means that D* 
correlates more strongly for that patient. 
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FIG. 2. Time course of two representative cases. The low- 
est trace shows the EEG and the second lowest the plasma 
concentration of the anesthetic estimated by measuring the 
concentration of the exhaled air. Depth of anesthesia is es- 
timated every 5 s by four different methods. The two top- 
most curves denote spectral parameters, the median frequency 
/median and the spectral edge frequency 95%, / ae 95 • The third 
curve shows the bispectral index (bis) as issued by the com- 
mercial monitoring equipment. The fourth curve gives the 
nonlinear correlation index D* . See text for discussion. 



thetic concentration between blood plasma (assumed to 
be equal to the concentration in the air at the end of ex- 
piration) and a hypothetical effect compartment (i.e. the 
brain or parts of it). It will be influenced by factors like 
brain perfusion, patients' age, amount of body-fat, etc. 
However, since the value of k is no longer needed once a 
monotonic relationship is assured, these influencing fac- 
tors need not be under close control for this approach. 

Based on the optimal values for k, for each case Spear- 
man's rank correlation coefficient p was computed (cf. 
Eq. ||). Table [0] shows that the highest correlation co- 
efficients were achieved for D* . In order to assess the 
statistical significance of this finding, we have performed 
a nonparametric test for the hypothesis that D* performs 
better than bis. Since the populations of values where ob- 
tained on the same patients, we evaluated the Wilcoxon 
matched-pairs signed-rank test (see Fig. ||) and found 
that the correlation is higher for D* with p < 0.05. 

For clinical applicability, it is not sufficient to have a 
good correlation for each patient individually. Rather, a 
monotonic relationship has to hold across patients. In or- 
der to illustrate the relevance of the relation between the 
brain concentration of Sevoflurane C e and its estimated 
EEG effect E for the whole group of patients, for each 
investigated EEG measure a scatterplot of C e against E, 
including results from all patients, is shown in Figs. |] 
and |5[ The corresponding rank correlation coefficients p 
are -0.06 for / mo dian, -0.33 for / SC 95, -0.63 for bis, and -0.72 
for D*. Due to the serial correlations within each trial, 
we don't attempt a formal test of the hypothesis that a 
(negative) correlation exists. Visual inspection as well as 
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FIG. 4. Cumulative plots of the correlation with estimated 
anesthesia depth of the two spectral indicators /median and 
/sc95, comprizing all considered cases. 

the value for p suggests that / mo dian is essentially useless 
as an indicator of depth of anesthesia in the investigated 
group of patients Also for the other spectral measure, 
/ S0 95, the correlation is rather unreliable. It is conceiv- 
able that in a study that excludes burst-suppression EEG 
or with additional preprocessing (e.g., by introducing a 
burst-suppression detector) the situation may be more 
favourable. Both bis and D* show a certain trend to 
decrease with increasing depth of anesthesia. Of all the 
studied measures, we find that D* comes closest to an 
overall monotonic relation to the depth of anesthesia. 

V. DISCUSSION 

In a relatively small but focused clinical study, we have 
compared four different methods to quantify depth of 
anesthesia by numerical analysis of EEG tracings. All 
four methods were defined outside the study whence the 
results can be considered to be out-of-sample. We have 
used the available data to answer the question which of 
the four measures comes closest to a monotonic relation- 
ship to the level of sedation. While for both of the two 




1 2 3 4 5 

C e [%] 

FIG. 5. Cumulative plots of the correlation with estimated 
anesthesia depth of the two nonlinear indicators, bis and D* , 
comprizing all considered cases. 

nonlinear measures, bis and D* , such a relationship seems 
to exist, the correlation is strongest for D*, despite the 
fact that bis had been specifically designed for this pur- 
pose in previous work. 

With the limited data base available, any interpreta- 
tion of the findings remains speculation. The success 
of D* in previous studies on the anticipation of epilep- 
tic seizures has been discussed in Ref. p(J , where it was 
pointed out that the connection with the epileptic alter- 
ations is lost if the temporal correlations are fully ex- 
cluded from the correlation sum, as one would have to 
do for a dimension estimation j|3]. It can be suggested 
that geometrical correlations are not essential but the 
ability to select structures - temporal or geometric - by 
amplitude, or length scale in phase space, makes the cor- 
relation sum superior to other nonlinear statistics. 

In earlier studies, D* has been found helpful for the 
anticipation of epileptic seizures |D| from intra-cranial 
EEG recordings. It has also been shown Jl6| that even 
on a cellular level D* could quantify the complexity of the 
synaptic input of a single neuron when inter-ictal epilep- 
tiform activity was induced in brain slices. The degree 
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of autonomous activity, and hence the number of degrees 
of freedom of a neuronal network, is thought to decrease 
during the development of epileptiform activity. Indeed, 
a reduction of D* extracted from intracellular recordings 
of single neurons was observed prior to visually detectable 
alterations ]l6| ]. 

In the common scenario of anesthesia, the number of 
degrees of freedom of human brain activity is lowered — 
as compared to the awake state — when nearly all neo- 
cortical neurons fall silent and only a central (thalamic) 
oscillator is able to evoke a cyclic excitation of the si- 
lenced neocortex ||. However, the number of degrees of 
freedom may still be huge and is probably not accessi- 
ble from an EEG time series It is even the more 
remarkable that D* correctly ranks the qualitatively dif- 
ferent EEG pattern of burst suppression as corresponding 
to the deepest level of anesthesia. At the onset of burst 
suppression, the EEG signature changes discontinuously 
with the estimated drug concentration, whence the rela- 
tively sharp transition in D* between the preceding deep 
sleep phases and burst suppression in Fig. ||. Still, from 
a biological point of view it remains unclear which elec- 
trophysiological properties contribute to the monotonic 
decline of D* under anesthesia — or epileptic activity. 

In summary, extraction of D* seems to be able to im- 
prove the quantification of depth of anesthesia from brain 
electrical activity, at least when Sevoflurane is used as 
anesthetic drug. For clinical applicability, other anes- 
thetics have to be investigated. The size of this study is 
too small to admit any optimisation of the algorithm. It 
is conceivable, however, that the clinical reliability can 
be enhanced, and the computational burden relieved, by 
a more thorough understanding of which structures are 
really picked up by the algorithm. 
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